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Abstract 

We consider the Falicov-Kimball model in two dimensions in the neutral case, i.e, 
the number of mobile electrons is equal to the number of ions. For rational densities 
between 1/3 and 2/5 we prove that the ground state is periodic if the strength of the 
attraction between the ions and electrons is large enough. The periodic ground state 
is given by taking the one dimensional periodic ground state found by Lemberger 
and then extending it into two dimensions in such a way that the configuration is 
constant along lines at a 45 degree angle to the lattice directions. 
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In the Falicov-Kimball model there are two types of particles. The electrons hop 
between nearest neighbor sites of the lattice while the ions do not hop at all. The electrons 
do not interact with each other except for the fact that they obey Fermi statistics. At 
most one ion is allowed at each lattice site and there is an on-site attraction between the 
ions and electrons. The Hamiltonian is 



<x,y> X 

where and Cx are creation and annihilation operators for the electrons. Vx is the 
occupation number for the ions, i.e., = 1 if there is a ion at x and Vx = if there is 
not. The sum over < x,y > is over nearest neighbor bonds in the lattice. (The factor of 4 
in front of the U is included for latter convenience.) The interaction of the electrons with 
the ions and the kinetic energy and fermi statistics of the electrons lead to an effective 
interaction for the ions. We will study the ground states of this model on the square lattice 
in the neutral case (the numbers of ions and electrons are equal) when the parameter U 
is large and positive. A review of rigorous work on the Falicov-Kimball model may be 
found in [Q. 

In any number of dimensions the ground state for density 1/2 is the checkerboard 
configuration for all f/ > In two dimensions with large U the ground states 

for densities 1/3,1/4 and 1/5 are known rigorously and are periodic 0, 0]. Numerical 
results indicate periodic ground states for numerous other values of the density ||15|. In 



one dimension Lemberger |]T0[ showed that for any rational density, the ground state is 
the periodic arrangement of the ions which is "most homogeneous." (He gives an explicit 
algorithm for determining the most homogeneous configuration.) His proof applies when 
U > Uq where Uq depends on the denominator of the rational density. 

The above results suggest that in two dimensions in the neutral case with large U the 
ground state is always periodic. However, Kennedy proved that for densities between 
1/5 and 1/4 other than 2/9 the ground state is not periodic. In this density range there 
is phase separation. The ground states for densites 1/5, 2/9 and 1/4 are periodic, but in 
between the ground state is formed by adjoining regions with these periodic configurations 
so as to produce the desired density. Haller proved that the same phenomenon occurs 
for densities between 1/6 and 2/11. So an analog of Lemberger 's result cannot hold in two 
dimensions for all densities. In this paper we prove that an analog of his result does hold 
for rational densities between 1/3 and 2/5. In this interval of densites the ground state 
is periodic. It is constant along lines of slope ±1 (one choice for the entire ground state), 
and its restriction to any horizontal or vertical line is the same as the one dimensional 
ground state found by Lemberger. We expect that this result is true for the entire density 
interval 1/3 to 1/2, but for technical reasons we cannot prove it for densities between 2/5 
and 1/2. 

There is a physical argument that suggests that the neutral model with rational density 
wants to spread out the ions as much as possible. Each electron spends most of its time at 
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Figure 1: The known results for the two dimensional neutral model with large U shown 
as a function of the density. The ground state is known to be periodic at all the densities 
at which there is a vertical line. In the interval (0, 1/6) the vertical lines indicate the first 
few densities l/iji^ + {n + 1)^) where the ground state is known to be periodic. 

a site with an ion since the attraction between electrons and ions is large. Now consider 
the effect of the kinetic energy of an electron. If nearby sites have ions, then the electrons 
at those sites will restrict the movement of the electron we are considering. So its kinetic 
energy is minimized by arranging the ions so as to maximize the space that each electron 
has to move in. However, Watson 0] emphasized that there will typically be a mismatch 
between the lattice and the ion configuration that maximizes this space in the absence 
of a lattice. Thus the lattice structure can frustrate the exclusion principle's attempts to 
put the ions in the "most homogeneous" configuration. In one dimension this frustration 
does not occur. In two dimensions Kennedy's and Haller's results on phase separation 
show it does occur for some intervals of density, while this paper shows it does not for 
other intervals. There does not appear to be any easy way to predict which intervals of 
density will have periodic ground states and which ones will have phase separation. 

The results for the neutral model on the square lattice with large U are summarized 
in figure |l|. The obvious open problem is to determine whether there is a periodic ground 
state or phase separation in these intervals where the behavior is unknown. However, the 
more important problem is to understand when there is phase separation vs. periodicity 
without the extensive calculations required in this paper or the papers on phase separation. 

The results on periodicity and phase separation in the square lattice depend on ex- 
tensive calculations of an effective Hamiltonian and thus depend heavily on the lattice 
structure. One possible way to gain some insight into the general periodic vs. phase 
separation question is to investigate other lattices. There are rigorous results on periodic 
ground states on the triangular lattice for some densities ^ . In three dimensions, nothing 
is known with the exception of density 1/2. 

The preceding discussion and indeed all of this paper is only concerned with ground 
states. For density 1/2 and any ?7 > it is known that there are at least two Gibbs states 
if the temperature is sufficiently low ^ and the number of dimensions is at least two. 
For densities 1/3, 1/4 and 1/5, a "quantum" Pirogov-Sinai theory must be used to prove 
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that for large U there are multiple Gibbs states at low temperature 0, |Tl|, |1^, |13|. The 
key ingredient in these approaches is a Peierls estimate. One can prove such an estimate 
for densities 1/6, 2/11 and 2/9 and so the methods should apply to these densities as well. 
For densities in (1/3, 2/5) we do not know how to prove the needed Peierls estimate, and 
so the existence of multiple Gibbs states at low temperatures is an open problem for these 
densities. 

We now turn to our main result. The precise statement is as follows. 



Theorem 1 Suppose p = p/q & [1/3,2/5] where p and q are relatively prime. We assume 
the lattice A is L x L where L is a multiple of 4q, and we impose periodic boundary 
conditions. Then there exists U{q) > such that for all U > U{q), the ground states 
are translations, reflections, and rotations of the configuration S defined by the following 
properties: S is constant along lines of slope +1, and the restriction of S to any vertical 
or horizontal line is the same as the one dimensional large U ground state with density 
P- 



First we give an overview of the proof. When U is large one may do perturbation 
theory in 1/U. We start by considering 3x3 blocks of sites and consider the perturbation 
series for the energy up to fourth order. Kennedy has shown that at this order, each block 
which looks like one of those shown in figure ^ contributes the same amount of energy, 
while any other block contributes a higher amount of energy. So if a configuration exists 
in which every 3x3 block looks like one of those from figure |^, it minimizes H through 
fourth order. It turns out that such configurations do exist. In fact, Watson WM has 
shown that they correspond to tilings by squares and parallelograms as shown in figure 
^. In such a tiling the configuration must be constant on lines of slope ±1, where the 
sign for the entire configuration is determined by the slope of the short sides of one of the 
parallelograms . 

Up to rotations and reflections, the ions in such a tiling must be one of the three 
types shown in figure ^. We may write the energy through eighth order in terms of the 
number of each vertex type. Through sixth order we find that every tiling has the same 
energy. At eighth order, tilings with only type B and C vertices have the same energy, 
but type A vertices increase the energy. So if there exists configurations that correspond 
to square-parallelogram tilings with no type A vertices, then they minimize H through 
eighth order. Such configurations do exist for densities in [1/3,2/5]. 

When we include the higher order terms, the most we can say at this point is that the 
ground state consists primarily of large regions corresponding to a square-parallelogram 
tiling with no type A vertices. Because of the lack of type A vertices, these regions must 
contain parallelograms. Furthermore, the slopes of the short sides of these parallelograms 
must be the same throughout each individual region. This means that within each region, 
the configuration must be constant on lines of slope -|-1 or —1, but not both, as it would 
in a region consisting entirely of type A vertices. 
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Figure 2: A typical example of a piece of one of the "+" regions. The thicker lines are 
the "stripes" of occupied sites. The dotted lines are inserted to show the squares and 
parallelograms. 

On each of these "+" and "— " regions, the problem essentially becomes one-dimensional, 
since we are now only concerned with the spacing between these "stripes" of occupied 
sites (see figure 0). It is tempting to try to apply a procedure to each region, similar to 
that of Lemberger for the one- dimensional case, in order to determine the cheapest stripe 
arrangement. However, there are two difficulties that arise. First there is the problem of 
the shape of these regions. Their boundaries may be very irregular whereas an extension 
of Lemberger's method would require applying periodic boundary conditions to a fairly 
regular shape such as a square or rectangle. Second there is a density constraint. The 
density inside the + and — regions need not be the same as the overall density. This 
is a problem because how large U needs to be in Lemberger's argument depends on the 
density. Even if we could say something about the configuration in each + and — region, 
this would not take into account the ions lying outside these regions, so we effectively lose 
track of the original ion density. We must therefore find more regular regions to work 
with that allow us to keep track of the original density. 

The first step in overcoming these difficulties is to partition the lattice into diamonds 
containing 8 sites, as shown in figure ^ If one of these diamonds intersects a 3 x 3 block 
not shown in figure ^, or it contains a type A ion, it may be associated with an increase 
in energy at eighth order. We call such diamonds "bad." The eighth order terms actually 
occur with f/^'', so the total energy cost of the bad diamonds is at least 0{U^'^)N, where 

is the number of bad diamonds. Each of the remaining "good" diamonds is contained 
in one of the + or — regions mentioned above. In a + region (— region), we group the 
diamonds into —45° strips (+45° strips) in such a way that the strips are bounded on 
each of the short sides by a bad diamond. This is shown in figure 0. 

If we were to apply periodic boundary conditions to a particular strip and compare the 
resulting energy with that of the energy of the strip in the original configuration, it would 
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seem that this difference would be proportional to the size of the boundary of the strip. 
However, consider one of the 45° strips of — diamonds. Inside the strip, the configuration 
is constant on lines of slope —1. These lines of constant configuration persist through the 
long sides of the strip, and continue in 5* until they encounter a bad diamond. So if there 
are no bad diamonds near the long sides, then imposing periodic boundary conditions on 
the strip is effectively the same as the boundary conditions the strip had to begin with. 
This allows us to show that the total difference between the energies of the strips in the 
original configuration and their energies with periodic boundary conditions is 0{U~^)N. 

We then glue all the strips of good diamonds together with the bad diamonds end to 
end to form one long thin strip. We apply periodic boundary conditions to the config- 
uration on this strip. The difference between the energy of this long strip with periodic 
boundary conditions and the energy of the original configuration is 0{U~^)N. Each bad 
diamond costs an energy 0{U~'^), so the energy of the original configuration is at least as 
large as that of our long strip. Finally, in the single long strip we can apply Lemberger's 
argument to conclude that the best arrangement is the most homogeneous one. 

Proof of Theorem |l]: We denote an ion configuration by S. We let H{S) be the 
minimum energy of H with this ion configuration, i.e., we minimize the Hamiltonian with 
respect to the electronic wave function. If U is sufficiently large, then the function H{S) 
may be expanded in powers of U~^. 

oo 

H{S) = U^'"'hra{S) (2) 

m=2 

(Details of this perturbation theory may be found in |l^.) Because we are on a square 
lattice, only the even values of m have a nonzero hm- hm{S) is a local function. It is 
a sum of terms, each of whose support is contained in the support of a closed nearest 
neighbor walk with m steps. We will let Hm{S) denote U^~"^hm{S). 

Let M be the number of 3 x 3 blocks in a configuration S which are not a reflection 
or rotation of one of those shown in figure Kennedy showed that 

aMU'^ < H2{S) + H^iS) + U'^L^a + bp) < ^MU'^ (3) 

where a > 0, 7 > 0, a, and h are constants. So after adding the constant U~^L^{a + bp) to 
the Hamiltonian, we see that every 3x3 block which is not shown in figure ^ contributes 
a positive amount of energy at fourth order (or lower), and the configurations in which 
every 3x3 block is one of those in figure |^ have zero energy through fourth order. 
The following lemma gives us a geometric interpretation of a configuration in which every 
3x3 block is one of those shown in figure ^. 

Lemma 1 (Watson [|14|| ) Let S be a configuration in which every 3x3 block is one of 
those shown in figure |^. Then there exists a tiling by squares and parallelograms as shown 
in figure^, in which every vertex corresponds to an ion in S. 
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Figure 3: Configurations which minimize the energy locally. 



Proof: Let S be a configuration satisfying the hypothesis of the lemma. Consider any 
occupied site. It is the center of a type 1, 2, or 3 block from figure ^. Now consider a 
3 by 3 block which shares six sites with the original 3 by 3 block we considered. Since 
the configuration at six of its sites is known and it must be one of the blocks in figure 
^, there are only a few possibilities for the configuration on the three unknown sites. 
By continuing in this manner and keeping track of all the cases one may determine the 
allowable local ion arrangements. These turn out to be rotations and refiections of the 
three arrangements in figure ^ g 




Type A Type B Type C 

Figure 4: Possible local ion arrangements for a configuration in which every 3x3 block 
is one of those in figure |^. 

From the lemma and equation (0), we see that the configurations which minimize 
the energy through fourth order correspond to tilings by squares and parallelograms. 
However, the densities one can obtain from such tilings is restricted. A configuration in 
which every ion is type C has density 1/3. One in which every ion is type A has density 
1/2. We may construct a configuration having any rational density between 1/3 and 1/2 
by placing an interface of slope ±1 between a region of type A ions and a region of type 
C ions, the size of each region being chosen to give the desired density. The ions on the 
interface will then be type B. So we see that from these square-parallelogram tilings, one 
may obtain any rational density in [1/3, 1/2]. (Of course, there are many other ways to 
do so besides the construction we have given.) 
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Consider a square- diamond tiling and let nA,nB,nc be the number of type A,B,C 
ions, respectively. These three integers do not determine the configuration, but they do 
determine enough of the local structure of the configuration that the energy H{S) may 
be computed through eighth order from uaiUb and rtc- We show in the appendix that 
for m = 6 and 8 

H^{S) = U^'^^iuAei + nse^ + nce^) (4) 

where 

= -3200, = 98000, = -1520, (5) 

ef = 33600, = 160, = -15120. (6) 

In a square diamond tiling, the numbers ua-, ub, and nc are not independent. Since the 
area is L^, we have 

5 

2nA + -riB + 3nc = L"^ (7) 
We are also keeping the number of ions fixed at pL^, so 

riA + riB + nc = pL^. (8) 
Using eqs. (|^ and (|) to solve for ub and ric in terms of ua, we have 

He{S) + Hs{S) - L'^f{U,p) = 15680f/-V. 

f{U,p) is a function of U and p, but there is no need to write it out explicitly. 

The preceding paragraph assumed that S was a tiling. In general, this will not be the 
case. Suppose S is not a tiling and M is the number of 3 x 3 blocks which are not rotations 
or refiections of those in figure ^ Then the number of sites in the region where S is not a 
tiling is at most 9M. So the terms intersecting this region contribute at most 0{U~^)M 
to Hq + Hs- Equations and (H) will no longer hold, but the difference between the 
right and left hand sides will be bounded by a constant times M. Thus for a general S, 
we have 

HeiS) + HsiS) - L'^fiU, p) = 15680f/- V + 0{U-^)M. (9) 

We may now determine the ground states through order 8. From inequality (^) we 
see that any configuration with M = will minimize the energy through fourth order. In 
equation (^, M appears at order 6. So for large enough U we want M = 0. Since the 
coefficient of ua is positive in equation (|]), we also want ua = 0. So the configurations 
which minimize the energy through order 8 are those where every ion is type B or C. These 
configurations correspond to square-parallelogram tilings with no adjacent ±45° strips of 
squares. They have zero energy through order 8 after we have subtracted L^/(f/, p). All 
we have shown so far is that zero is a lower bound for the energy through eighth order for 
densities in [1/3,2/5]. We can construct a configuration that attains this lower bound as 
follows. A region with only type B vertices has density 2/5, while a region with only type 
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C has density 1/3. By adjoining a type B region with a type C region with an interface 
between them at 45° we can obtain any density in [2/5, 1/3] and attain our lower bound. 

If we just consider the energy through eighth order, then any configugation that has a 
"bad" 3x3 block or a type A ion will have energy higher than the ground state energy. 
The crucial point is that there exists a positive constant c so that difference of the energy 
and the ground state energy is at least cf/~^ times the number of these "local defects" in 
the configuration. We now move on to the higher order terms. 

Recall that p = p/q & [1/3, 2/5] where p and q are relatively prime, and the lattice A 
is L X L where L is a multiple of 4g. We partition A into diamonds containing 8 sites as 
shown in figure |[ 



Figure 5: Partitioning of the lattice into diamonds containing 8 sites. 

We define a diamond h to be had if it contains a type A ion or intersects a bad 3x3 
block. Otherwise, h is good. We let N denote the number of bad diamonds. A good 
diamond is one which is contained in a region of the configuration corresponding to a 
square-parallelogram tiling with no type A ions. We know from the above that the energy 
through eighth order is bounded below by a term proportional to the number of bad 3x3 
blocks plus the number of type A ions, which is proportional to N . So inequality (^) and 
equation (|^) tell us that 

H2 + H^ + Hq + Hs + U-^L\a + hp) - L^f{U, p) > cNU-' (10) 

where c > is a constant. For future convenience, we define 

H = H + U-^L'^{a + hp) - L'^fiU, p) 
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Let S he a configuration with density p{S) G [|,|]- We label a parallelogram in S 
with a + if its short sides have slope +1, and a — if its short sides have slope —1. Now, 
every good diamond must contain a type B or C ion, both of which are the vertex of a 
parallelogram. Any other parallelogram intersecting this diamond must have the same 
sign. Thus the good diamonds may be labelled + or — in a natural way. Furthermore, 
since all the ions in a good diamond are type B or C, two diamonds that share an edge must 
have the same sign. One may think of 5* as being divided into regions where the diamonds 
are +, regions where the diamonds are — , and regions of bad diamonds. A + region and 
a — region must be separated by a region of bad diamonds. Inside every + region, S is 
constant on lines of slope +1, and inside every — region it is constant on lines of slope 
— 1. We call this the stripe property. More, however, may be said about the arrangement 
of the occupied stripes inside each good region. Since every ion must be type B or C, 
each pair of occupied stripes must be separated by one or two unoccupied stripes. Also, 
since there are no type A ions, we may not have consecutive stripes corresponding to the 
following occupation sequence: occupied, unoccupied, occupied, unoccupied, occupied. 
This rules out the checkerboard configuration, the sign of which is ambiguous. When a 
region has the stripe property with no type A ions, we say that the stripes are correctly 
spaced. 

We now divide up the — regions into 45° strips as follows. Choose any good diamond 
whose sign is — . As shown in figure ^ we extend this set in the directions ±(1, 1) until we 
hit a bad diamond. Every diamond in this 45° strip is good, and, since one of them is — , 
they must all be — . We may then carry our labelling one step further and label this strip 
with a — . Note that every good — diamond is contained in exactly one such — strip. 




Figure 6: A good "— " diamond is chosen. The sites in the shaded region are contained 
in good "— " diamonds. 

We divide up the + region into —45° strips by changing all the signs in this construction. 
We label these strips with a +. See figure |^ 
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bad diamonds 



Figure 7: Left: a — strip. Right: a + strip. The configuration inside the — strip (+ strip) 
is constant along fines of slope —1 (+1). 



We now take tfie strips and rotate tfiem so tfiey are all fiorizontal and tfien glue tfiem 
togetfier to form one long strip. Now rotate tfie bad diamonds and glue tfiem onto one 
end to form an even longer strip. Tfiis strip will contain all tfie sites of tfie original A. 
Tfiis is sfiown in figure Tfie strips fiave tfie stripe property. After rotation tfiis means 
tfiat in figure ^ tfie portion of tfie configuration tfiat came from tfie strips is constant 
along vertical lines of lattice sites. In figure ^ vertical lines of lattice sites only contain 
two sites, so tfiis simply means tfiat eitfier botfi of tfiese sites are occupied or botfi are 
empty. In tfie portion of figure ^ wfiicfi contains tfie bad diamonds, tfie configuration may 
not satisfy tfiis property. However, in tfie statement of tfie tfieorem L was assumed to be 
a multiple of 4g. Tfius tfie total number of ions is a multiple of 4p. Tfiere are an even 
number of ions in eacfi + and — strip, so tfiere must be an even number of ions in tfie bad 
region of our tfiin strip. Tfierefore we may rearrange tfie configuration in tfiis portion of 
tfie strip so it also fias eitfier zero or two ions in eacfi vertical line of lattice sites. We now 
glue tfie sfiort sides togetfier and tfie long sides togetfier to form a torus, cafi it T. Let 
denote tfie resulting configuration on T. We want to bound \H{S) — H{Si)\. As we sfiall 
see, we will only need to estimate tfiis difference at orders > 10. 




bad 
diamonds 



Figure 8: Tfie tfiin strip obtained by gluing togetfier tfie + and — strips as well as tfie 
bad diamonds. 

Tfie perturbation tfieory leads naturally to an expression for if as a sum over nearest 
neigfibor closed walks 7. 

7 
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Lemberger gave a nice formula for the contribution of a walk: 



where I7I is the number of steps in 7 and m(S', 7) is the number of ions at sites in the 
walk 7. Obviously, Hn is equal to the sum over walks with n steps. 

We want to set up a correspondence between the terms in the perturbation series of 
H{S) and H{Si). So we need to set up a correspondence between walks on the torus T 
and the original square A. We think of Z'^ as the covering space of the torus T. Then a 
walk 71 on T has a natural lift to a walk 71 on Z^. Each site in T comes from a site in 
A. We let 7 be the translate of 71 such that 7(0) is the site in A that corresponds to the 
site 7i(0) in T. The map 71 — > 7 takes walks on T to walks on A. 

There is a problem with this correspondence. The lifted walk 71 need not return to 
its starting point. In fact the lifted walk will return to its starting point if and only if 
the original walk 71 is homotopic to zero. So we redefine H{Si) and H{S) to be the sum 
of all the terms in the perturbation series that come from walks homotopic to zero. If a 
walk on A is not homotopic to zero, then it has at least L steps. (L is the length of the 
square A.) Thus this changes H{S) by terms of order U~^. So the total change in H{S) 
is of order L'^U^^. This is tiny even compared to the energy cost of a single bad diamond. 
The change in H{Si) is not so small since the torus T is quite narrow in one direction. 
But we are free to define H{Si) any way we wish. It merely serves to organize our proof. 

Now suppose that 71 — > 7 and 7 does not encounter a bad diamond in S. Then 7 lies 
entirely in either a region of + diamonds or a region of — diamonds. Consider the first 
case. So the strips have slope —1 and the configuration is constant along hues with slope 
+ 1 in the vicinity of the walk 7. It follows that 7 sees the same configuration in S that 
71 sees in Si, and so uj{'y) = Thus the terms in Hn{S) — Hn{Si) that do not cancel 

come from walks 7 in 5* that do not encounter a bad diamond. So the number of terms 
that do not cancel is of order n^N. Hence we have 

\Hr,{S) - Hn{Si)\ < U^-^d'n^N (12) 

Note that although satisfies the stripe property, its stripes may not be correctly 
spaced. The idea now is to form a finite sequence of configurations, {Si\f^-^, such that 

K — 0{N), and the stripes in the final configuration 5";^: are correctly spaced. The 
strategy is as follows. The reason why the stripes in 5*1 may not be correctly spaced is 
because of the introduction of the bad diamonds, as well as the cuts. There arc 0{N) 
stripes in the bad region of Si and 0{N) cuts, so we may hope to correct the spacings 
with 0{N) modifications, each modification costing at most C/^~"c" at order n. 
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In making these modifications we must be sure that each one fixes a "fault" in Si 
without creating a similar fault. All the modifications follow the same basic procedure. 
We illustrate this by considering a pair of adjacent occupied stripes in Si. Since the 
density of 5*1 is < 2/5, if there is a pair of adjacent occupied stripes, there must be a 
pair of adjacent unoccupied stripes. As shown in figure ^ we introduce two "markers," 
one between the pair of adjacent occupied stipes and one between the pair of adjacent 
unoccupied stripes. We let W be the portion of the configuration between the markers. 
With the two axes of rotation a and b as shown in figure ^, we form by rotating W 
about the axis a. We then replace W by and, if necessary, rotate about the 
axis b so it fits into place. (Whether or not we need to rotate about b depends on the 
number of sites in W.) S2 is then the modified configuration. Any walk which gives a non 
zero contribution to Hn{S2) — i/„(S'i) must cross one of the markers, and so must visit 
one of the 8 sites adjacent to the markers. So we have that \Hn{S2) — Hn{Si)\ < [/^""c". 
Our modification has allowed us to separate a neighboring pair of occupied stripes 
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Figure 9: Separating a pair of neighboring occupied stripes. W is rotated about axis a 
and, if necessary, axis b. 

without creating a new pair. So we see that we can get rid of all such pairs with 0{N) 
modifications. 

The next problem we focus on is the possibility of more than two consecutive unoccu- 
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pied stripes. If this occurs, since the density is > 1/3 there must be at least one pair of 
occupied stripes separated by a single unoccupied stripe. Figure |10| shows how we may 
reduce the number of adjacent unoccupied stripes by one by a procedure similar to the 
above. 
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Figure 10: Reducing the number of consecutive unoccupied stripes. 

After 0{N) repetitions of the above procedures we reach a configuration in which 
every pair of occupied stripes is separated by one or two unoccupied stripes. We now face 
the possibility that our configuration contains type A ions. These ions could have come 
from 5*1 or they could have been formed by the modifications above. But we only made 
0{N) modifications, and 5*1 contained at most 0{N) type A ions. So our configuration 
can contain at most 0{N) type A ions. We let a "1" denote an occupied stripe and a "0" 
denote an unoccupied stripe. The stripe sequence signifying a type A ion is then: 1, 0, 1, 
0, 1. If this occurs, then since the density is < 2/5 there must exist the following stripe 
sequence: 1, 0, 0, 1, 0, 0, 1. In figure |Tl| we see how the procedure used above can get rid 
of these type A ions. At first glance, it may seem that this last modification, although 
removing a type A ion at the left end of W in Si, may have created a type A ion at the 
right end of in S'j+i. However, if we do get a type A ion in S'j+i after the inversion of 
W, it must have been there before the inversion, in Si. Thus the number of type A ions 
is indeed reduced. 

The last configuration in this sequence of modifications, Sk, has correctly spaced 
stripes. Furthermore, K = 0{N) and |if„(S'j+i) — Hn{Si)\ < U^~'^c^ for each consecutive 
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Figure 11: Removing the type A ions. 



pair in our sequence. This gives us the following telescoping sum 

K-l 

H{S) - H{Sk) = [H{S) - H{S,)] + - if (13) 

1=1 

The construction of our sequence above and inequality (|^) imply that 

J2 \HniS) - Hn{SK)\ < K U'-^c^n^ < cNU-'' (14) 

even n >10 even n >10 

In the last inequality we used the fact that K is of order A^. Since Sk satisfies the stripe 
property and its stripes are correctly spaced, we know that 

8 

Hn{SK)=0 

even n=2 

Equation (|T3D and the bounds ([1^) and ([M]) imply 

H{S) > H{Sk) + aNU-^ - cNU'\ (15) 

Consider the torus T on which Sk is defined. Let Sl be a configuration on T which 
satisfies the stripe property and such that the restriction of the configuration to a line in 
one of the lattice directions is the one- dimensional ground state for density p{S) when U is 
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large as determined by Lemberger. {Sl is unique up to lattice symmetries.) We claim that 
among all the configuration on T having density p{S) and satisfying the stripe property, 
the energy is minimized by Sl- In particular, the claim implies H{Sk) > H{Sl)- 

To establish the claim we show that there is a close relationship between the pertur- 
bation series for the two-dimensional model and that of the one-dimensional model. The 
terms in the Hamiltonian at order 2m may be written in terms of the walks as 

H2m= Yl ^(^) (16) 
7: |7|=2m 

with ^(7) given by eq. (]TT|). For the one- dimensional model the only change that needs 
to be made in this formula is that 7 is only summed over walks in the one-dimensional 
lattice. We will show that if the configuration has the stripe property then 

7:|7|=2m,,7 hor 

where 7 hor means that the sum is over walks which only take steps in the horizontal 
direction. So these walks are one-dimensional although they may start at any point in 
the two-dimensional lattice. The combinatorial factor c(m) is given by 

k=l ^ ^ 

To show (|T^ we think of a walk as a sequence of 2m letters chosen from L (left), R 
(right), U (up) and D (down). The number of L's must equal the number of R's, and the 
number of U's must equal the number of D's. (This is because of our redefinition of the 
Hamiltonian so that the perturbation series only contains terms coming from walks that 
are homotopic to zero.) Let Pu) be the projection of the walk that is obtained by replacing 
each U by an R and each D by an L. (This assumes that we are in the case where the 
slopes are —1. For slope +1, we replace U by L and D by R.) Then for configurations that 
satisfy the stripe property, H{u) = H{Puj). Now let u' be a walk with only L and R steps. 
Since the walk is closed there are m R's and m L's. We obtain all walks u with Pu = uj' 
by replacing some number of L's with D's and the same number of R's with U's. If k is the 
number of L's that we replace, then there are ('^) ways to make the L — D replacements 
and an equal number of ways to make the R — U replacements. So the number of uo with 
Puj = uj' is c(m) as defined above. Equation (0) says that for stripe configurations the 
two-dimensional Hamiltonian is the same as the one-dimensional Hamiltonian order by 
order except for a multiplicative factor at each order. It is straightforward to check that 
Lemberger's argument applies unchanged when the c(m) factors are included. 

The stripe arrangement in follows the Farey tree, just as the ion arrangement in 
the 1-d case. So Sl consists of a unit cell of length g, which is periodically extended to 
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cover the torus. Our original square lattice A has sides whose lengths are multiples of 4g. 
So we may build a configuration on A which has the same stripe arrangement as Sl, 



and satisfies H{Sy,) = H{Sl)- We may then rewrite inequality (|T5[) as 



H{S)> H{St) + aNU~'^ -cNU-^. (19) 

Assuming U > \fcfa, this shows that 5*^ is a ground state. To see that these are the only 
ground states, we consider two cases. If iV > then the above bound says li{S^ > H{Sf;). 
If = 0, then S satisfies the stripe property with the stripes correctly spaced. So we 
know from Lemberger's argument that H{S) > H{S^), with equality holding only when 
5* = 5*^ up to a lattice symmetry. | 
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Appendix 



In this appendix we derive eq. @ which expresses the sixth and eighth order terms in 
the Hamiltonian in terms of the number of ions of types A,B and C. Recall that the 
Hamiltonian at a given order is given in terms of walks by 



H2m = 5Z 



7:|7|=2m, 



This may be rewritten as 



H2m — Cm,xS{X) 



(20) 



(21) 



X 



where the sum is over subsets of the lattice and S{X) is 1 if the configuration S has an 
ion at every site in X and otherwise. 

Recall that I{S) denotes the sets of sites at which there is an ion for the configuration 
S. If S{X) is to give a non zero contribution to Hq{S) + Hg^S), X must be a subset of 
I{S) which is contained in a closed nearest neighbor walk of six or eight steps. Let S 
correspond to a tiling by squares and parallelograms. Then the only subsets X C I{S) 
which can contribute at orders 6 and 8 are those shown, up to rotations and reflections, 
in figure |12|. We will refer to an ion and the eight ions arranged about it as in figure 
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Figure 12: All possible subsets of I{S) (up to lattice symmetries) which contribute to the 
energy at orders 6 and 8 when S corresponds to a square-parallelogram tiling. The lines 
serve to distinguish the subsets. 
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^ as the "local arrangement" of the ion. Recall that the local arrangement about an ion 
in I{S) must look like one of the three shown in figure So for any X which gives a 
nonzero contribution to Hq{S) +Hs{S), there is at least one ion x such that X is a subset 
of the local arrangement about x. There may be more than one such x. For example, if X 
consists of a single ion, X will be contained in 9 local arrangements. To compensate for 
this over counting, we define Cm,x = '^(x) where i^{X) is the number of ions whose local 
arrangement contains X. Now fix a type A ion and let L denote the local arrangement 
about it. Then for even m < 8 we have 

~ ^ ] Cm,X 

This may be written as 

xefigure 13; 

XQL '-' 

where "X G figure means that X is a rotation or refiection of one of the arrange- 
ments from figure and ta,x is the number of translates of X appearing in the local 
configuration about a type A ion. and are given by the same formula with L taken 
to be the local arrangement of a type B or C ion respectively. 

Table |l] gives the coeffecients Cq^x and Cg^x for each X in figure 0. The last three 
columns give the numbers ^fxj' ^^"^ example, that eg = 

64 X 1 - 384 X 2 + 96 X 2 - 674 X 4 = -3200? 

We should note that one can use the above type of argument to show that an equation 
of the form (|D must hold. The values of e^, and e';!^ can then be computed by carrying 
out the perturbation theory for Hm{S) for several different periodic configurations S and 
then using the results to set up a linear system of equations for e^, and e^. We have 
done this as an independent check on the calculation above. 
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A 


C6,X 


C8,X 


A 

A 


D 


L- 


1 


64 


112 


1 


1 


1 


2 


-384 


768 


2 


3/2 


1 


3 


96 


-1360 


2 


1 





4 


216 


-768 





1 


2 


5 


24 


448 





1 


2 


6 


-672 


11520 


4 


2 





7 




512 


4 


1 





8 




32 


2 








9 




192 


2 








10 




1152 


2 


1 


1 


11 




768 


2 


1 


1 


12 




-1728 





2 


4 


13 




-5184 





1 


2 


14 




-384 


8 


2 





15 




-1248 


4 


1 





16 




3840 


4 


1 





17 




23040 


1 


1/2 





18 




3840 


4 


1 
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960 


1 









Table 1: Values of Cm,x used to calculate the energy of types A, B, and C ions at orders 
6 and 8. 
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